f 

A TWO-DIRECTIONAL TARGET OPTIMIMTION MODEL 



V 



Gregory R* Hamel in 



SS5SSS”“ 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

A TWO-DIRECTIONAL TARGET OPTIMIZATION MODEL 

by 

Gregory R. Hamelin 
March 1979 

Thesis Advisor: Gilbert T. Howard 

Approved for public release; distribution unlimited. 



118863 ^ 




i 



I 



i 



I 

I 



rUDLEVKNOXUUriARV 

F I ED 



SeCUniTY CLASSIFICATION OF THIS ^ACC O^mn 



Dmtm Bnffmd) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. NCPONT NUMSCR 


2. OOVT ACCESSION NO. 


2. NCCIPlCNT*S CATALOG NUMSCN 


4. riTLt (and SuMlIm) 

A Two-Directional Target Optimization 
Model 


5. TYNE OF NEPONT 4 PEPlOO COVEPEO 

Master's Thesis; 

March 1979 


S. PCNFONMING ONG. NCPONT NUMSCN 


7. AUTHONfa; 

Gregory R. Hamelin 


S. CONTRACT OR GRANT NUMBCRfai 


S. FCNFONMINO ONOANIZATION NAMC ANO AOONCSS 

Naval Postgraduate School 
Monterey, California 93 94 0 


iO. program CLCMCnT, PROJECT, TASK 
AREA A WORK UNIT NUMBERS 


11. CONTNOLLING OFFICC NAMC ANO AOONCSS 

Naval Postgraduate School 
Monterey, California 93940 


12. PEPOPT OATE 

March 1979 


11. NUMEEP of PACES 

61 


14. mONITONINO agency name 4 A00NESS<l< dl/f«rant tram Contnlllng Otilea) 

Naval Postgraduate School 
Monterey, California 93940 


18. SECURITY CLASS, (at thia ripen) 

Unclassified 


\%m. OECLASSIFICATION/DOWNGRAOING 
SCHEDULE 



16. distribution STATCIACNT rol thim Rmport) 



Approved for public release; distribution unlimited. 



17. OISTNtSUTION STATCMCNT (oi thm mhmtrmct mtimrmd In MSoak 30, U dUItmti from Mmport) 



IS. su^plemcntany notes 



It. KEY WONOS (Continum on eow^mm &Idm It nncmmmmrr Idmntlir br Sloe* numkor) 



20. ASST pact (Ccniimim mt rormmn 9idm It nnenmmmr and Idmniltr tr tioeir numkor) 

This paper presents an algorithm for computing the optimal 
target path for two aircraft traversing a target area from 
different directions. There are constraints on the maneuver- 
3-t>ility of each aircraft which prohibit it from attacking 
every target. The algorithm chooses a subset of targets whose 
iestruction will yield maximum value to the attacking force. 



DO 1473 COITION OF 1 NOV SS IS OSSOLCTC 

(Page 1) 0102-014- 6601 I 



TTTJrT.A.q.qTPT-FD 

ItCUUlTY CLASSIFICATION OF THIS FAOt (Whan DM tnttaW) 







* 




UNCLASSIFIED 

tacu»wTv eLA»ii»ie*Tio>i ow this n««« 



The basis of the algorithm is the branch and bound method, 
with upper bounds computed by dynamic programming. Several 
variations are considered, such as payload limit, an increased 
number of aircraft from each direction, and a three-directional 
attack. An, example problem is solved using the basic model. 

A Fortran IV computer program is included. Computation 
time versus problem characteristics is discussed. 



DD Fornj 1473 
S/?} *{fl‘?)2-ni4-6601 



UNCLASSIFIED 



2 



euAUiricATioN 



I 



I 



I 

I 

I 

I 



Approved for public release; distribution unlimited. 



A TWO-DIRECTIONAL TARGET OPTIMIZATION MODEL 



by 



Gregory R. Hamel in 
Lieutenant, United States Navy 
B.S., United States Naval Academy, 1972 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN OPERATIONS RESEARCH 



from the 

NAVAL POSTGRADUATE SCHOOL 
March 1979 



^ I 

c ( 



ABSTRACT 



This paper presents an algorithm for computing the optimal 
target path for two aircraft traversing a target area from 
different directions. There are constraints on the maneuver- 
ability of each aircraft which prohibit it from attacking every 
target. The algorithm chooses a subset of targets whose de- 
struction will yield maximvim value to the attacking force. 

The basis of the algorithm is the branch and bound method, 
with upper bounds computed by dynamic programming. Several 
variations are considered, such as payload limit, an increased 
number of aircraft from each direction, and a three-directional 
attack. An example problem is solved using the basic model. 

A Fortran IV computer program is included. Computation 
time versus problem characteristics is discussed. 
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I . INTRODUCTION 



Optimal assignment of targets to aircraft on a strike 
mission can greatly increase the effectiveness of an air attack 
against ground installations. The purpose of this paper is to 
present an algorithm which will, given the location and rela- 
tive military/industrial worth of key enemy positions, select 
a sequenced subset of targets whose destruction will yield 
maximum value to the attacking force. It assumes that the 
decision maker has complete knowledge of the targets and their 
value. A stochastic extension can easily be incorporated 
when considering hardened targets. This is done by multiply- 
ing the probability of killing a particular target by its value. 

The initial coiorse of an aircraft approaching a target 
area is denoted by (j>. Due to anti-air defenses, there is a 
limit placed on the maximum number of degrees which an aircraft 
may deviate from its initial course. This angular deviation is 
denoted by 6. Its effect is to restrict the aircraft's movement 
at any point to a cone whose vertex angle, 20, is bisected by 
the aircraft's initial course, 4>. Thus at any point (x,y) in 
the target area, the aircraft's course options are between 
4i-e and <I)+6. This creates an ordering among the targets which 
dictates the sequence in which they must be considered. The 
allowable deviation is pictured in Figure 1. 
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Not Possible 




Not Possible 



(x,y) 



Figure 1. Cone of Allowable Course Deviations 

Because of' the course constraint, there are only a limited 
number of paths through the target area which are feasible. 

The task of the decision maker is to choose that set of tar- 
gets which yields the optimal combined value. 

References [1] and [2] present a method for determining 
the optimal set of targets for one aircraft traversing an area 
containing MN targets. The problem is formulated as an MN+1 
stage dynamic programming problem [Ref. 3]. The targets are 
numbered in decreasing order from an imaginary line drawn tan- 
gent to the boundary at which the aircraft enters the target area. 
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The stages are represented by straight lines drawn through 
the targets/ parallel to the boundary tangent line. Stage n 
corresponds to the parallel line drawn through target n. 

The starting point of the aircraft, outside the target area, 
is stage MN+1. It is located so that any target is accessible 
from stage MN+1. The stage diagram is depicted in Figure 2. 
Should two targets be equidistant from the boundary tangent 
line, either one may be assigned the next sequential number, 
and the problem reduces to MN stages . 
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Figure 2. The Stage Diagram for One Aircraft 

Entering from below the Target Area . 
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It is important to realize that were the aircraft to enter 
the target area, say, from the right, the stage diagram would 
appear as in Figure 3. 



STAGE 




Figure 3. The Stage Diagram for One Aircraft Entering Target 
Area from the Right. 
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The state variable denotes the lateral position of the 
aircraft at stage n. is the decision as to the heading of 

the aircraft as it moves from stage n to stage n-1. is 

restricted to lie in the set of feasible headings S^, from 
(j)-6 to (()+0. The state variable is then a function of 
and This function is referred to as the stage trans- 

formation t. 

The return function for stage n is denoted by r^. Letting 
be the lateral position of target n, and be its value. 
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The problem is then written as 



Maximize 


MN 

Z r (X ) 

n=l ^ ^ 




Subject to: 


^n=^(^n+l'°n+l^ 


ri“l f • • • f 






n=2 f • • • ,MN+1 



An efficient algorithm for solving this problem is dis- 
cussed in Section VI, but the problem can also be easily solved 
graphically. To do so, the course deviation angle must be 
viewed from the perspective of the stage diagram. Figure 4 
shows whether target n is feasible for various values of 
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n+1 



^n+1 

not feasible 




feasible 



^n+1 

not feasible 



Figure 4. Cone of Feasibility 

Using Figure 2 as an example, the maximum return possible 
for stages (n-1) to 1 is recorded on each stage line n, 
n = 1, MN+1, for every possible value of x^^. Assuming 

the value of each target is one, Figure 5 gives the solution 
for one aircraft entering from below the target area. 

Tracing back from stage MN+1, the optimal set of targets 
is 8, 7, 6, 4, 3, 2, for a value of six. 

Dynamic programming can extend this problem to M aircraft 
by increasing the number of state variables to M [Refs. 1, 2, 
and 3] , as long as all M aircraft are attacking from the same 
direction. But when the aircraft enter the target area from 
different directions, a pure dynamic programming solution is 
no longer possible, since the association between the target 
and the stage is no longer valid. A target in stage three 
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for one aircraft may be in stage ten for another, and the re- 
cursive equations have no meaning. 




Figure 5. Graphical Solution to the One-Aircraft Problem. 
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This paper presents a solution to the problem of two air- 
craft attacking a target area from different directions. A 
model is formulated, and the algorithm developed. An example 
problem is included. Also included is a computer program for 
implementing the algorithm and a discussion of its effective- 
ness. Because the terminology becomes quite involved, a 
glossary of terms is provided in Appendix A. 
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II. PROBLEM FORMULATION' 



In formulating this problem, it is assumed that if an 
aircraft attacks a particular target, that target is destroyed. 

In other words, there are no misses. As was stated in Section 
I, this does not preclude a stochastic approach for hardened 
targets where each has a probability of being destroyed. 

For ease of notation, the two aircraft are denoted air- 
craft A and aircraft B. 

The targets are niimbered as in Section I, but now each 
target has two niombers, one with respect to aircraft A, and one 
with respect to aircraft B. They are denoted A^ and B^ , respec- 
tively, i= 1, ..., MN, j = 1, ..., MN. It is critical to real- 
ize that if any target A^^ has the same coordinates as any tar- 
get Bj, then that A^ is the identical target Bj . In fact, for 
every there will be a Bj identical to it. Figures 2 and 3, 
which are identical target areas, should clarify this point. 

The value of A^ will be denoted by V[A^], and the value of 
Bj by V[Bj] . 

A feasible set of targets for a single aircraft is one in 
which the course required to go to each successive target is 
within the allowable course deviation cone described in Sec- 
tion I. PA (I) denotes one of the feasible sets of targets 
for aircraft A. It is convenient to think of PA(I) as a 
path of targets which aircraft A will attack. Any PA (I) com- 
pletely disregards aircraft B in that it is computed as if only 
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one aircraft were attacking. Similarly, PB(J) denotes one 
of the feasible sets of targets for aircraft B. 

The value of a set of targets, or a path, is the summa- 
tion of the individual target values comprising that set. The 
values of PA(I) and PB(J) are denoted by V[PA(I)] and 
V[PB(J) ] , respectively. 

In order for the solution of a two-aircraft problem to be 
feasible, the sets of targets for aircraft A and aircraft B 
must be mutually exclusive, that is, they must have no targets 
in common. This is so because one aircraft is sufficient to 
destroy the target. No additional return would be realized 
by the other aircraft attacking the same target. If this 
feasibility constraint were not required for optimality, aircraft 
A (or aircraft B) might forego the opportunity to attack other 
targets in order to attack the target that both aircraft have 
in common. 

The optimal solution is achieved when PA (I) and PB(J) are 
chosen so as to 

Maximize V[PA(I)] + V[PB(J)] 

Subject to: PA(I) D PB(J) = J0. 
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•III. THE ALGORITHM 



The branch and bound method forms the basis of the solution 
algorithm for the two-aircraft problem. It is discussed in 
theory in Refs. 4, 5, 6 , and 7. Only a small fraction of the 
possible solutions to the problem is actually enumerated. The 
remaining solutions are eliminated from consideration through 
the application of bounds that establish that such solutions 
cannot be optimal . 

The algorithm begins by considering all possible combina- 
tions of paths for both aircraft. It then breaks this set of 
all possible combinations into smaller and smaller subsets, and 
calculates for each an upper bound on the value of the best 
paths contained therein. The bounds determine the partitioning 
of the subsets and eventually identify an optimal path for both 
aircraft. The branch and bound method represents the subsets 
as nodes of a tree and the partitioning of the subsets as a 
branching of the tree. 

Node one consists of all possible combinations of paths for 
both aircraft. Using the single aircraft dynamic programming 
(D.P.) method of Section I, the optimal path for aircraft A 
is computed for the direction from which A enters the target 
area. This path is denoted by PA(1) . By the same method, the 
optimal path for aircraft B is computed for the direction from 
which B enters. This resultant path is denoted by PB(1). It 
is important to realize that both PA(1) and PB(1) are computed 
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as single aircraft optimizations, and that all MN targets are 
possible elements of PB(1) , even those in PA(1) . 

If PA(1) and PB(1) have no targets in coinmon, then the 
paths form the optimal solution to the two-aircraft problem, the 
value of which is V[PA(1)] + V[PB(1)], and the algorithm stops. 
However, if there are targets in common, then as was mentioned 
in Section II, the solution is not feasible. The summation 
V[PA(1)] + V[PB(1)] represents instead an upper bound, denoted 
by UB(1) , for the optimal solution. The summation of the 
values of the points of intersection of PA(1) and PB(1) is de- 
noted by VINT(l) . The target in the intersection which has 
the highest value is denoted by X^(l) . This target could be 
written in terms of aircraft A or aircraft B. When the dis*» 
tinction is necessary, X^(l) will be written as either ^^(1)^ 
or X,„(l)_, for aircraft A and aircraft B, respectively. Geo- 
graphically, however, X^(l), X^(l)^, and X^^CD^are identical. 

The goal is to find optimal paths for aircraft A and B 
which have no targets in common. The nature of the D.P. solution 
for a single aircraft path is such that it seeks out those fea- 
sible targets with the highest value. In node one, both A and 
B sought X^g(l). If the set of path combinations in node one 
was restricted so that aircraft A had to take X^(l)^ and air- 
craft B could not take X^(l)g, the same target, then that 
point of intersection would be eliminated. But perhaps the 
optimal solution requires B, not A, to take X-_(l). It could 
even require that neither path include X^g(l). Therefore, 
to include all possibilities, the set of node one is broken 
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into two subsets, one which requires that the path of aircraft 
A exclude X^(l) , giving aircraft B the option of taking it 
or not, depending on the single aircraft D.p. solution for B. 
The other subset requires that the path of aircraft B exclude 
X,„(l), giving aircraft A the option. Thus node one branches 
to form nodes two and three. 

The restrictions placed on each aircraft at node I are de- 
noted by R(I) . If, for example, R(I) = A., B., B,, the path 

1 3 K 

of aircraft A would be required to exclude target A^, and the 
path of aircraft B would be required to exclude targets Bj 
and Bj^. 

The branching, with restrictions, is illustrated in 
Figure 6. 




AB 



(1)b 



Figure 6. The Start of a Tree 

Now consider node two with the restriction vector R(2) 
containing the single element X^(l)^. Again using the 
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single aircraft D.P. method, only this time with the restric- 
tion that the path of aircraft A not include target X^(l)^, 
the optimal path for aircraft A is computed. The restriction 
can be incorporated into the solution techniques of Section I 
by temporarily assigning a large negative value to X^g(l), 
thereby making it highly unattractive as an element of the 
optimal path. Denote the resultant path PA (2) , and return 
the original value to X^g(l) • As was noted, restricting 
PA (2) to exclude X^(l)^ places no restrictions on aircraft 
B. Therefore, PB(2) will be identical to PB(1). Just as 
with node one, UB(2) , VINT(2), and ^^(2) can be computed. 

Node three is considered next, with restriction vector 
R(3) containing the single element Since this means 

that the path of aircraft B must exclude X^(l)g, a large 
negative number is assigned to X^(l)g,PB(3) is computed, 
and the original value is returned to X^(l)g. Since R{3) 
places no restrictions on PA(3) , it is identical to PA(1) . 
Next, UB(3), VINT (3), and ^^(3) are computed. 

A terminal node is one from which branching may still 
occur. Nodes two and three are terminal nodes. Since only 
two branches may emanate from any node, node one is no longer 
terminal, and need not be considered further. 

The next step in the algorithm is to choose the terminal 
node which has the highest upper bound. Assume it is node 
J (at this point, J is either two or three, whichever has the 
higher upper bound) . Node J then branches to form nodes four 
and five. The restriction vector R(4) will equal R(J) with 
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The restriction vector 



the addition of element 
R(5) will equal R(J) with the addition of element X._(J)_. 

Again the restricted paths are computed and the algorithm 
continues until the stopping condition is met. 

The branching is always done in pairs, a left branch and 
a right branch, as was shown in Figure 6. No valid upper 
bound comparisons can be made until both branchings have 
been performed and the upper bounds of both new nodes have 
been computed. 

The stopping condition for the algorithm occurs when, 
following a double branching, a terminal node is found whose 
upper bound is greater than or equal to all other terminal 
node upper bounds, and whose paths for aircraft A and aircraft 
B are mutually exclusive, that is, they have no targets in 
common. This solution is optimal because its value is equal 
to the upper bound of that node, making it the best solution 
possible for that node. And since its value is at least as 
good as the best solution possible for all other terminal 
nodes, it is a global optimum. 

The steps of the algorithm are summarized in Figure 7 
at the end of this section. 

There is one other calculation which can be made at each 
node which is of some interest. The lower bound for each 
node I is denoted by LB (I) and is equal toUB(I) minus VINT (I). 
This is a lower bound because if all the points of intersection 
on PA(I) and PB(I) were given to either aircraft A or air- 
craft B, the value of the resultant solution would be LB (I). 
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If any terminal node I has a lower bound on the optimal value 
which is greater than or equal to the upper bound of any other 
terminal node J, then node J may be completely dropped from 
consideration. It is said to be fathomed, and is no longer 
terminal. This will not speed up the algorithm or reduce the 
number of branchings required, since the algorithm would 
never branch from node J anyway. However, if computer storage 
space were critical, it would be advantageous to incorporate 
lower bounding, since once a node was fathomed it could be 
removed from storage. It will not be employed in this al- 
gorithm. 

The algorithm guarantees that an optimal solution will be 
found. However, it suffers from a limitation common to all 
branch and bound methods. For any untried problem, it is im- 
possible to tell beforehand exactly how much computation will 
actually be necessary to find the optimal solution. Depending 
on the way the problem is set up, it could converge to the 
optimum very quickly, or for a large, difficult problem it 
could require such excessive branching that it becomes compu- 
tationally prohibitive. This would be the case if the allow- 
able course deviation were very large, the angle between at- 
tackers very small, and there were multiple optimal paths. 
However, in sample problems of target optimization, the al- 
gorithm converges very quickly, as will be shown in Section 
VI. 
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